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Abstract 

Classical motion of positronium embedded in a magnetic field is studied, and computation shows 
the emergence of chaotic orbits. Recent work investigating quantum behavior of this system pre- 
dicts extremely long lifetimes @. Chaos assisted tunneling however may lead to significant 
shortening of the lifetime of this system. 
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Two interacting charged particles placed in a magnetic field exhibit chaotic motion. This 
has been studied for the hydrogen and Rydberg atom [|I| and the scattering of electrons on 
positive nuclei [Q.This has an impact on the electrical conductivity of fully ionized plasmas 

i i- 

Here we study the classical motion of positronium in a magnetic field. In the absence of 
a magnetic field the positronium has a very short lifetime. It was found by Ackermann et. 
al.[^ that in a strong magnetic field the positronium can have an extremely long lifetime 
"up to the order of one year" 

We find that the classical motion is chaotic, which usually leads to chaos assisted tunneling 
which should significantly reduce the lifetime of this system. 

The calculation includes the case of crossed electric and magnetic fields, provided that 
the ratio of the field strengths E/B does not exceed the speed of light. In this case the 
electric field can be eliminated by a Lorentz transformation. 

The motion of two particles with charges +e and -e of equal mass m moving, in a 
uniform magnetic field B are described by the equations 

. ri-r2 
mri = evi x is (11 

Atxeq \r1-r2y 

mr2 = -era x B + -r (2) 

47reo \r1-r2 

Adding (|l]) and (||) gives 

m(fi + f's) = e(ri - 7-2) x B (3) 
Introducing the new coordinates 

ri + r2 = R (4) 
ri - 7-2 = r (5) 

and integrating Eq. results in 

rriR — er X B = a (6) 
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where ck is a constant vector. Subtracting Eq.(|]) from ([l|) gives 



mr = eR x B (7) 



and using Eq.(||) 
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e c T 
mr = — (er x B + a) x B (8) 

Introducing the cyclotron frequency Uc = eB/m, and choosing B pointing in the z 
direction B = Be^, Eqs. (§) and (H) become 

f/ujl = (r X 63) X 63 + a X 63/(6^) - ^ (9) 



it/uJc — r X e^ = ct/eB (10) 

With the dimensionless variables Uct — > t and r {271 e^B"^ /niY^^ r, one arrives to the 
dimensionless equations of motion 

f = (r X 63) X 63 + a' X 63 — r/r^ (11) 

R-r X e3 = a (12) 

where a = a/eB is the dimensionless constant vector. Since r x has no component 
in the z direction, it is convenient to chose a coordinate system where the initial value of 
Rz{0) = 0, so the constant cx is a vector in the x — y plane, ex = aei + 662. Without loss 
of generality one may chose either a=0 or b=0. So Eq.([TTD becomes 



r + r± = —r/r^ + bei (13) 
where _L means perpendicular to the z axis. 



First we study the two dimensional case 



X + X = —x/r^ + b 

y + y = -y/r^ 



(14) 
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with the Hamiltonian 



H = + - 1/r - bx (15) 

Since energy is conserved phase space is three dimensionaL The potential energy V = 
— 1/r — bx is singular at r — 0, and develops a minimum when b > 1.89. |l] shows 
equipotential lines for b=3. Surface of section plots in the x, px plane have been computed for 
different values of the two parameters: b and the energy E. p.a^ hows four plots, E = 0,b = 1 
and b = 3] E = —.5, 6=1 and 3. The plots appear regular without any chaotic orbits. 
We have carried out many more computations for different values of E and b with similar 
results. It appears therefore that an additional constant of motion exists, but the analytic 
expression has not been found. We have also carried out computations to find the largest 
Lyapunov exponent which turned out to be zero as expected for non-chaotic orbits. 

Turning to the three dimensional case we study the a = limit, in this case it is 
convenient to introduce polar coordinates where the system is described by the Hamiltonian 

H = P;/2 + Pl/2 + Pl/{2p') + pV2 - 1/ V7 + ^ (16) 

where = a;^ + , and ip is an ignorable coordinate so Pip = const. This gives the 
equations of motion 

P 

z 

(17) 

dp — P^ ^ipl P (^2+^2)3/2 

dH z 

dz (p2+^2)3/2 

A surface of section plot (p, Pp) in the z=0 plane is shown in ^ for E=-.5, = .25. The 
existence of chaotic orbits is obvious, so the three dimensional equations of motion are not 
integrable. To show that chaotic orbits exist for o: 7^ 0, the largest Lyapunov exponent has 
been computed for the three dimensional case for b=3, as shown in ^, using the algorithm 
as described in Ref. It converges to a value larger then zero as expected. 



dH 
dPp 

dH 
OP, 



p. = - 



4 



In conclusion the motion of positronium immersed in a magnetic field is chaotic in the 
classical limit, therefore the long lifetime predicted in the quantum limit is unlikely. 
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FIG. 2. a: Surface of section plots in the x — Px plane,E=0,b=l 
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FIG. 2.b: E =0,b = 3 
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FIG. 2.d: E = -0.5, b = 3 
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FIG. 3: Surface of sections plot for the three dimensional case, where b= in the p — Pp plane 
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FIG. 4: Computation of the largest Lyapunov exponent, E=-0.5, b=0.1 with initial conditions 

XI = 0.65,yi = 0.0,zi = 0.0 

P,i = O.Q,Pyx = 0.25,P,i = 1.31222067 

X2 = 0.65000001,^2 = 0.0,^2 = 0.0 

P^2 = 0.0,Py2 = 0.25,P^2 = 1.31222064 



